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Static three quark potential in the quenched lattice QCD 
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We study the static three quark potential using lattice QCD simulation. At zero temperature, the three quark 
potential is extracted from the baryonic Wilson loop and fitted to the phenomenological form, the sum of the 
Coulomb term, linearly rising term and the constant. We compare two types of the linear term, "Y" and "A" 
type, and find that the string tension almost coincide with uqq in the former case. We also investigate the three 
quark potential at finite temperature using the Polyakov loop correlation. In the deconfined phase, the screened 
potential is observed similarly to the case of QQ system. 
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1. Introduction 

Determination of the static quark-antiquark 
(Q-Q) potential has been one of most vital sub- 
jects of the lattice QCD simulation On the 
other hand, the status of the three quark (3CJ) po- 
tential is much less settled, in spite of its impor- 
tance in the baryon spectroscopy. In this paper, 
we present our recent result of lattice calculation 
of the static three quark potential . 

It has been widely used the static Q-Q po- 
tential expressed as the sum of the short range 
Coulomb term and long range confining linear 
term. In the case of the 3Q potential, two dis- 
tinct treatments of the linear term are frequently 
applied. In this paper they are denoted as "A- 
type" and "Y-type" forms. The former expresses 
the linear potential as the sum of two-body po- 
tentials, which are proportional to the distance 
between each pair of quarks. The latter treats 
the potential with the string linking three quarks 
at the point with which the total length of the 
string takes the minimum length, L min . Y-type 
ansatz has been used in the phenomenological 
studies with the string tension a^Q ~ a QQ> an( i 
successfully describes the gross feature of the 
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baryon spectrum [Q. A-type ansatz is, however, 
also used in some model calculations with the 
string tension a a effectively around half of &qq ■ 
The lattice simulations have been done mainly 
in 1980's ^-?J, and there is one recent prelimi- 
nary work [B|. Among them, several authors in- 
terpreted their results as the support of the A- 
type picture ■ This embarrassing situation 

is mainly caused by the insufficient extraction of 
the ground state contribution from the baryonic 
Wilson loop. In nineties, calculation of Q-Q po- 
tential became much reliable by making use of the 
smearing technique which enhances the ground 
state contribution to the Wilson loop This 
procedure is directly applicable to the three quark 
system, and is expected to enable us reliable ex- 
traction of the ground state contribution. 

Our first goal of this work is an extraction of re- 
liable static 3Q potential using the smearing tech- 
nique. With the lattice result, we discuss prefer- 
able description of the potential. These subjects 
are described in the next section. 

Another important subject on the 3Q potential 
is the thermal effect on it at finite temperature. 
Again the lattice computation is much retarded 
compared with the Q-Q potential. We present 
our preliminary result at T > T c in Section 3. 
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Figure 1. Baryonic Wilson loop. 

2. Zero temperature result 

The three quark potential is extracted from the 
baryonic Wilson loop, 
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W 3Q (f u r 2 ,r 3 ,T) = ^e abc t a , b , d Ufufuf, (1) 



Pexp{ig / dx^A^x)}, 



(2) 



where P denotes the path-ordered product along 
the path i\. (fc =1, 2, 3) as shown in Figure |l|. 
The position of the contraction of three lines Uk 
can be chosen arbitrarily, and the result should 
not depend on this point. We choose the con- 
tracting point as the spatial origin, and locate 
three static quarks on each three spatial axis. The 
same constellation is averaged making use of the 
translational invariance. 

The numerical simulation is done on the lattice 
of size 12 3 x 24, with the standard Wilson gauge 
action at (3 — 5.7 in the quenched approximation. 
The lattice spacing is determined as a ~ 0.19 fm 
by setting the string tension of Q-Q potential to 
be a = 0.89 GeV/fm. The (baryonic) Wilson loop 
is measured on 210 configurations with applying 
the smearing technique to enhance the ground 
state contribution. The potentials are extracted 
from the (baryonic) Wilson loop by the single ex- 
ponential fit in appropriate fit regions, where the 
Wilson loop is dominated by the ground state 
contribution 

Let us start with the Q-Q potential. It is widely 
accepted that the lattice result of the static Q-Q 
potential is well described by the form 

V(r) =C + or. (3) 

r 

In Table |, we list the results of the fit in three 
kinds of fit ranges: on-axis fit, fits with off-axis 



data in the whole r region and at r > 2. Since 
our lattice spacing is rather coarse, the fit with 
the off-axis data may suffer from the rotational 
symmetry breaking effect. In this sense, the fit 
with the on-axis data is preferable to determine 
the parameters. For the 3Q potential, however, 
we inevitably treat the off-axis constellations, and 
then it would be meaningful to examine the typi- 
cal size of the rotational symmetry breaking effect 
in the Q-Q system. For this purpose, we carry 
out the fit including the off-axis data with two fit 
ranges. The fit using all data (all r in Table re- 
sults in large x 2 value, and this is mainly caused 
by the result at the short distance region where 
the statistical error is smallest and the cutoff ef- 
fect is most severe. We alternatively fit at r > 2, 
which is more appropriate for comparison with 
the result of the 3Q potential. In this case the fit 
results in x 2 l^df = 3.3, and we regard this value 
as a guide of the rotational symmetry breaking 
effect on this lattice. In both cases, the fit re- 
sult of the Coulomb coefficient much differs from 
the value of on-axis fit than the case of the string 
tension. 

Figure || shows the lattice result of the 3Q po- 
tential from the viewpoint of the y-type ansatz. 
The horizontal axis, L m i n , is the minimum length 
of the flux lines linking three quarks. Three flux 
connecting the junction with quarks take the an- 
gles 27r/3, which balances three forces with equal 
strength, in the case that all three angles of the 
"quark triangle" are less than 27r/3 [|J. The data 
are fitted to the form 

V{ri,f 2 ,r 3 ) = C Y -A Y (— + — + — 

III I II • (4) 

The result is listed in Table |l[ This form of the 
potential well describes the observed result, and 
the string tension is almost same as in the case 
of Q-Q potential. The x 2 of the fit is in the same 
order as in Q-Q with fit range r > 2. Considering 
the rotational symmetry violation effect at this (3, 
this seems to be a consistent result. 

Alternative ansatz, A- type, is represented with 
the fitting function 



F(ri,f 2 ,f 3 ) 



C A -A. 



A 



1 



1 

^23 



r 31. 



3 



C A a X Z /N df 


3Q(Y) 


0.914(20) 0.1316(62) 0.1528(20) 


4.0 


3Q(A) 


0.934(20) 0.1405(60) 0.0858(15) 


11. 


Q-Q 
on-axis 
all r 
r > 2 


0.629(16) 0.2793(12) 0.1629(47) 
0.686( 7) 0.3351(55) 0.1545(20) 
0.696(24) 0.395 (35) 0.1561(38) 


0.59 
16. 
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Table 1 

Fit results of the Q-Q and the 3Q potentials. 



+o-a(^i2 + r 2 3 +r 3 i), 



(5) 



and the obtained parameters are also listed in 
Table [l]. As the result of fit, we find that the 
string tension is much smaller than the Y-type 
fit, although the constant and the coefficient of 
the Coulomb term are almost the same as the Y- 
type results. In the fit, a a is slightly larger than 
the half of ay. Geometrically, 2L m i n ~ L\ = 
T\i + + roughly holds in most cases. This 
explains that cta — \ a QQ- There is another pos- 
sible physical explanation for a a — \ a QQ using 
the SU(3) C Casimir scaling. For the Coulomb 
coefficient, one gluon exchange brings the color 
factor 1/2 in the perturbation theory as the ra- 
tio of the SU(3)c Casimir factor. The linear 
term is however nonperturbative and then there is 
no clear explanation to keep the SU(3) C Casimir 
scaling beyond the perturbation. In addition, the 
Y 2 of the A-type fit is much larger than the case 
of the Y-type fit, where the goodness of the fit is 
comparable with Q-Q potential case with r > 2. 
Although the present simulation may suffer from 
rather large finite a effect, our result shows that 
the Y-type ansatz is more preferable with the lat- 
tice data. 

Our conclusion at this stage is the following. 
Y-type ansatz is more preferable for the lat- 
tice result although A-type ansatz is not com- 
pletely excluded with present systematic uncer- 
tainty. Adopting Y-type picture, we have found 
a universal feature of the string tension as a 3 Q ~ 
uqq. To confirm this result and investigate the 
3Q potential in more detail, simulations at larger 
/3's are in progress. Preliminary result on the 
16 3 x 32 lattice at j3 = 6.0 is consistent with the 




Figure 2. 3Q potential from the viewpoint of 
Y^-type ansatz. L m i n is the minimum length of 
the flux linking three quarks. 

result at (3 = 5.7. 

3. Finite temperature result 

Now we turn on the temperature and observe 
how the potential changes with the thermal effect 
|J. In this paper, we focus on the static 3Q po- 
tential in the deconfined phase. Above the crit- 
ical temperature T c , the quarks and the gluons 
are expected to be liberated and form a plasma 
state, in which the potential is expected to be 
screened. Then one expects that the screened po- 
tential takes the Yukawa potential like as the Q-Q 
potential fl(i|| , and genuine three-body potential 
(in Y-type picture) would disappear. 

On the finite temperature lattice, the static po- 
tential is extracted from the Polyakov loop corre- 
lation, 



P3(ri,r 2 ,r 3 ) = (P (n) P (r 2 ) P (r 3 )) 
= c ■ exp(-Y 3 Q(ri,f 2 ,r 3 )^T), 



(6) 



where the Polyakov loop P{x) is defined as 

P(x) = Trn^o" 1 u ± The Q-Q potential 

is also measured in similar way. 

To simulate the finite temperature system, we 
adopt the anisotropic lattice |]llj. An advantage 
of the anisotropic lattice for our study is that it 
enables us to work out several values of temper- 
ature without varying the lattice cutoff. As the 
gauge field action, we use the Symanzik improved 
action at the tree level |T^] . 

The numerical calculation is done on the lat- 
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tices of sizes 16 2 x 24 x N t at [3 — 4.56 and 
7g = 3.45, in the quenched approximation. At 
Nt = 96, which is realized as the zero tem- 
perature lattice, we determine the renormalized 
anisotropy as £ = 3.95(2) from the Wilson loops 
in the spatial-spatial and temporal-spatial planes. 
The static Q-Q potential at T = gives the cut- 



offs 



= 1.61(1) GeV and a" 1 = 6.36(5) GeV. 



The critical temperature is estimated by the rapid 
change of the Polyakov loop and its susceptibil- 
ity, and is found between N t = 24 and 25, and 
slightly near to N t — 24. For the investigation 
of the static potentials in the deconfmed phase, 
we generate 60 configurations at Nt = 20, which 
roughly corresponds to T ~ 1.22T C . 

Figure || shows the Q-Q potential and the 3Q 
potential in the deconfmed phase. At the long 
distance, The Polyakov loop correlators approach 
to the multiples of the expectation value of single 
Polyakov loop: P 2 -> (P){P)* and P 3 -> (P) 3 . 
These contributions are subtracted from the po- 
tentials. The Q-Q potential above T c is known 
to behave as the Yukawa- type potential (l(J. In 
present calculation, the Debye screening mass of 
the Q-Q potential strongly depends on the fit 
range and is estimated around 0.5-1 GeV. At 
this Nt, we have measured the 3Q potential with 
only equilateral constellations. In Figure ||, we 
show the 3Q potential divided by 3, and find good 
agreement with the Q-Q potential, 

^Q(ri2=r 23 = r 31 =r)~3V/ Q0 (r) (T>T C ). (7) 

Above T c , since the quarks are deconfmed, there 
is no reason to assume that the combination of 
quarks are in the color-singlet states. The ques- 
tion is whether the 3Q potential can be repre- 
sented as the sum of two body potential. Our 
present result shows that the Vsq/3 behaves as 
almost same as the Q-Q potential at least for 
the equilateral 3Q constellations. For more defi- 
nite conclusion, we need precise results including 
other 3Q constellations and a systematic compar- 
ison between VJjq and Vqq at various T. 

The calculations have been done on NEC SX-4 
at Research Center for Nuclear Physics, Osaka 
University and NEC HSP at INSAM (Institute 
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Figure 3. Q-Q (open symbols) and 3Q (filled) 
potentials at N t = 20 (T > T c ). The 3Q potential 
is divided by 3. 

for Nonlinear Science and Applied Mathematics), 
Hiroshima University. 
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